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Abstract 

The studied model describes a particle that obeys a one-dimensional non- 
linear Schrodinger equation in the potential of a double-well. Transitions 
between the two lowest self-trapped states of this system under the influence 
of the external time-dependent perturbation are studied in the two-mode ap- 
proximation. If the perturbation dependence on time is harmonic with the 
frequency oj, then transitions between the states become possible if the am- 
plitude of the perturbation F exceeds some threshold value F c (oS); above the 
threshold motion of the system becomes chaotic. If the perturbation is a 
broadband noise, then transitions between the states are possible at arbitrar- 
ily small F and occur in the process of the system's energy diffusion. 
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I. INTRODUCTION 



The nonlinear Schrodinger equation (NLS), 



dip 



A?p + U{r)^ + \\^\ 2 ilj, 



(1) 



2m 



serves as a tool for the account of many physical phenomena. The stationary version of this 
equation was introduced by Deigen and Pekar to describe the self-trapped (autolocalized) 
states of electrons in deformable crystal lattice ( |J, see also @). Gross || and Pitayevskii [3| 
derived Eq. (1) as the mean- field approximation for the macroscopic wavefunction if) (r, t) of 
the Bose - Einstein condensate of the non-ideal Bose gas at vanishing temperature (see also 
H). The third avatar of the NLS came in the realm of the nonlinear optics, where ip (r,t) 
represents the envelope of a quasimonochromatic electromagnetic wave. The exact solution 
of the one-dimensional homogeneous (U (f) = 0) NLS found by Zakharov and Shabat || 
formed the modern paradigm of the soliton theory 0]. 

Although the solutions of the NLS were extensively studied, it seems that comparatively 
little is known about their properties in case of time- dependent potentials. In this paper we 
address ourselves to a specific problem of this class. We shall consider the one-dimensional 
Eq. (1) with the potential that consists of two parts: a permanent potential U (x) that has 
the form of symmetric double well, and some time-dependent potential V (x, t) that will be 
called perturbation. In the absence of perturbation the properties of the stationary solutions 
of Eq. (1), that have the form 



depend essentially on the nonlinearity parameter A. At small A there is an infinite set of 
modes (2) that have symmetric wave functions - odd or even, $ (x) = ±$ (— x). With the 
increase of A at some threshold value of A c a pair of stationary solutions (2) with broken 
symmetry, $ s (x) ^ ±$ s (— x), comes into being. These solutions describe the states of the 
particle that are self-trapped in one of the wells of the permanent potential. For A < the 
corresponding energy E s is lower than that of any symmetric mode |§. 




(2) 
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Our main concern will be the following question: if the initial state of the system is one 
of these self-trapped states, then how the system will evolve under the influence of the non- 
stationary perturbation? In particular, can the perturbation transfer the system completely 
to the opposite self-trapped state? 

There is a favorable circumstance that allows to simplify the problem. It happens that 
at moderate A > A c even the self-trapped modes of high asymmetry can be accurately 
described by linear combination of the two lowest symmetric modes, the even $o (x) and 
the odd $1 (x) |§ . Therefore studying the problem we can restrict ourselves by analysis of 
the evolution of the two-level system. In Section II we derive the basic equations for this 
model. In Section III we study the influence of the external perturbation that harmonically 
depends on time. The evolution of the system under the influence of the broadband noise 
is studied in Section IV. Section V contains the concluding discussion. 

II. THE BASIC EQUATIONS 

For the future use we introduce the following quantities related to the (supposedly real) 
eigenfunctions $o (x) an d 3>i (x): 



Let's represent the wave function of the system by the superposition of two lowest symmetric 
modes, 



where $ = h" 1 + XJu). By substitution of Eq. (4) in Eq. (1), consequent multiplication 
by $i (x) , and integration over the coordinate x we get the following system of two equations 
for the complex amplitudes bf. 





V> (x, t) = b (t) % {x) e" w + h (t) $i (x) e~^\ 



(4) 
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(5) 



ih^ = -\b \b \ 2 J 00 - 2\b \h\ 2 J Q1 - Xb* o bp O i^°- 0l)t 
at 

+b V 00 (t)+hV 01 (t) e^-ft)*, 

rib 

ih-± = -A&i \h\ 2 J u - 2A6x |6 | 2 Joi - Afc 2 ,^ Joie^ "^* 

+hVu (t) + b V Q1 (t) e -^-^\ 
where the matrix elements of the perturbation are given by the integrals 

Vij (t) = f $i (x) V (x, t) (x) dx. (6) 



The system (5) conserves the norm of the state ip (x, t) (the sum of probabilities \bo\ + \bi\ = 
1), and the common phase of the wave function (4) is physically irrelevant. Therefore we can 
describe the evolution of the system by just two real variables. The complex amplitudes could 
be cast in the form 6j = ^fnl exp(— i^), where rij and $j are real time-dependent variables. 
Let's introduce the population difference A = n — n\ and the phase 6 = 2 (i? — $i) + 
2 (/3 — A) £. Then the system (5) turns into equations 

A = -B(l - A 2 ) sin 6 + F (t) y/1 ~ A ' sin -, 

2 2 (7) 

A 

= -Q + 2AA + 2BA cos - F (t) , cos YG (t) , 

VI - A 2 2 w ' 

where the following notations have been introduced: f2 = 2 (/3i — /?o) + A/iT 1 (Jn — Joo), 

A = A/T 1 (4J i - Joo - Jn)/2, 5 = A^Voi, F(t) = 47rV i (t) and G (*) = 

2h~ 1 (Voo (£) — Vii (t)). All these quantities have the same dimensionality namely that of the 

frequency. If the nonlinearity parameter vanishes, A = 0, then Eqs. (7) become equivalent 

to the well-known Bloch equations [jJJ. 



The nonlinear Bloch equations (7) can be considered as a pair of canonical equations for 
the conjugated variables A, of the non- autonomous system with one degree of freedom 
with the Hamiltonian function H = H + Hi (t), where the unperturbed Hamiltonian is 

H = ttA + AA 2 - B(l- A 2 ) cos©, (8) 

and the perturbation is 




(9) 



In what follows we refer to the value of the function H as the energy of the system. In 
the absence of perturbation the system (7) has two trivial stationary solutions, A = ±1 
and arbitrary, that correspond to the symmetric eigenstates <3>o and $i respectively, and 
two non-trivial fixed points A = — 0/2 (A + B), = or = 2ir, that correspond to the 
pair of self-trapped states that below will be called the stationary states. These states are 
divided from the bulk of the phase space by a separatrix (see Fig. 1). 

For the future numerical calculations we need to specify the parameters of the unper- 
turbed Hamiltonian (8). We have chosen the following set of values: Q = 5.388, A = 1.902, 
and B = 2.022. With this choice the stationary states which are located at A = —0.686 
correspond to the minimal energy of the system E_ = —3.871, the separatrix coincides with 
the isoenergetic line E — E s — —3.486, and the maximal energy E + = 7.290 corresponds to 
the line A = 1. 

III. THE HARMONIC PERTURBATION 

We shall assume that the even (diagonal) perturbation is absent, G(t) = 0, and the odd 
(non-diagonal) has the form 



The numerical experiments show, that for given values of frequency u and initial phase of 
the perturbation there is a threshold value F c (cu, 0) of its amplitude such that for F < F c 
the phase trajectory of the system remains indefinitely within one loop of the separatrix, 
whereas for F > F c the phase trajectory crosses the separatrix, does it repeatedly and may 
come close to the opposite stable point. The dependence of F c (cu, 0) on the initial phase is 
feeble: the relative variations of the thresholds due to the change of the initial phase have 
the order of a few percent. Hence for the time being we ignore this dependence and shall 
speak only about the dependence F c (cu) (see Fig. 2). 



F(t) = F sin (ut + 0) . 



(10) 



5 



The abrupt change of the character of motion at a certain threshold value of the per- 
turbation magnitude strongly indicates on the onset of the global stochasticity that comes 
from the overlap of resonances [JIT] , |I^| and the destruction of the noble tori |13|,|T2j. This is 
indeed the case: by taking some phase point at the separatrix for the initial conditions, one 
can see that at the same threshold values F c (u) the stochastic layer around the separatrix 
explodes and covers the vicinity of the stable states. However, even at rather small excesses 
of F over the threshold the crossing of separatrix comes fast, after about ten periods of the 
field. At these times the chaotic nature of the system's dynamics remains concealed: the 
motion seems regular and rather simple. Therefore we can try to explain the behaviour of 
the dependence F c (u) in the frame of regular dynamics. 

Specific features of the perturbing terms in Eqs. (7) create technical complications 
that are irrelevant to the nature of the phenomenon. The main qualitative features of the 
separatrix crossing under the influence of the harmonic field could be explained with a toy 
model - the one-dimensional Duffing oscillator with the equation of motion 

x + x — x 3 = F sin ujt (11) 

and the initial conditions x (0) = 0, x (0) = 0. This model also has a stable point, surrounded 
by a separatrix. 

We assume the frequency detuning 5 = to — 1 to be small, \5\ « 1. If the perturbation 
F is weak, then the non-linearity of the oscillator could be neglected, at least in the lowest 
approximation. Then the solution of Eq. (11) has the approximate form 



F . 5 
x [t) — - sin —t cos 
o 2 



5 ' 

ut--t 



(12) 



By assumption that that the oscillator could be linearized in all the range \x\ < 1, from this 
law of motion we find a crude estimate for the threshold of the separatrix crossing, namely 
F c =|*|. 

Now we improve this estimate by taking into account the nonlinearity. Let's represent 
the motion of the oscillator in the form x (t) = A cos (ut + <p), where A and (p are slowly 
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varying functions of time. The law (12) corresponds to the equations of motion for the slow 
amplitude and phase, 

F 5 
A = --cosy?, V? = ~2> ( 13 ) 

with the initial conditions A (0) = 0, ip (0) = ir. Let's now replace the eigenfrequency of 
the oscillator f2 = 1 in the RHS of the second of Eqs. (13) by the eigenfrequency of the 
nonlinear oscillator Q (A) that depends on the amplitude. For the model (11) for small A 
we have Q (A) = 1 — 3A 2 /16. Consequently the evolution of the system could be described 
by the system of equations 

A = -^cos<p, ( p = J--^-A 2 , (14) 

with the same initial conditions. The threshold of the separatrix crossing could be found 
from the condition that oscillations may reach the saddle points: maxA(t) = 1. From the 
second of Eqs. (14) it is seen, that if 5 > 0, then the phase shift decreases monotonously, thus 
decreasing the rate of the amplitude growth. If S < — 3/16 = —0.188, then the phase shift 
increases monotonously while the amplitude stays below its critical value, A < 1; and again 
the rate of the amplitude growth decreases with time. But in the band —3/16 = —0.188 < 
5 < the phase shift at first grows, then reaches the maximum and starts to decrease, 
passing the zero value at some later time to- Consequently, there are two moments in which 
the amplitude growth rate is maximal, t = and t = t . Thus one may expect that the 
dependence F c (u) will have the minimum somewhere in the range 13/16 = 0.812 < uj < 1. 
The numerical solution of Eqs. (14) shows that this is true: the minimum of F c (uj) is reached 
for lu = 0.87, about the middle of this band (see Fig. 3). 

To find the condition for the separatrix crossing, the equations (14) should be solved 
on the finite interval of time, while the phase shift reaches the value ±ir/2. This could be 
done in different ways, that may produce analytical estimates for the threshold values. We 
restrict ourselves by the only example. In the exact resonance (at 5 = 0) in the zeroth 
approximation the amplitude dependence on time is linear, A = Ft/2. By substitution of 
this expression in the RHS of the second of Eqs. (14), we have in the first approximation 
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tt/2 

/ r 2/3 cos ede 



0.0666 (16) 



Mt) = K-^F 2 t\ (15) 

The time t m when the amplitude A(t) reaches the maximum could be found from the con- 
dition ifi (t m ) = tt/2. Hence from the first of Eqs. (14) in the first approximation we have 
the threshold value of the perturbation in the exact resonance: 

= 

v ; 16 

o 

It agrees with the result of the numerical solution of the system (14) with the accuracy about 
6%, but differs from the value obtained in the numerical experiment by a factor about 1.5. 

The studied model (11) shares with the system (7) the "soft" character of the nonlinearity 
of oscillations around the stable points: the eigenfrequency of oscillations decrease with the 
growth of their amplitude. This common feature is responsible for the similar behaviour of 
F c (u) in two models (compare Figs. 2 and 3) - namely, the presence of non-zero minimum 
at a frequency somewhat lower than that of the small oscillations, CIq. 

Now we return to the case of the chaotic motion of the system above the threshold. For 
the system with the Hamiltonian H = H (A, 0) + V (A, 0) sinu;t with a small perturbation 
V the energy half-width of the stochastic layer is given by the Melnikov - Arnold integral 



CO 



AE = I l^A+^e)smutdt, (17) 



— oo 



where A (t) and (t) are taken for the unperturbed motion on the separatrix [11] . Thus 
we can expect that the motion in the phase space will persist inside the domain limited by 
the isoenergetic line E = E s + AE. Since for the Hamiltonian systems the phase volume 
is conserved, we may expect the invariant density in the phase space to be uniform. This 
leads to the invariant distribution of the energy values w (E) of the form 

w(E) = r] nTE)> (E^<E<E S ), 

1 ) (18) 
W {E) = V n(E)> (E S <E<E S + AE), 
where rj is the normalization constant, and the factor "2" in the first line accounts for the 
double degeneracy of the energy states. The comparison of this distribution with the one 
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obtained in the numerical experiments is shown in Fig. 4. The general agreement is clearly 
present, in spite of rather large value of the perturbation magnitude. The discrepancy 
between the distributions for the energy values around and above E s + AE is due to the 
borderline resonances of the stochastic layer and could be anticipated. 

If we define the vicinity of the stable state by the condition E < E*, then the average 
fraction of time spent in this domain is 

f dE (E,-E) . . 

and the average transition time from vicinity of one of the self-trapped states to vicinity of 
the other is about 

(20) 



m(#*) 

where r is the energy relaxation time. For the over-threshold perturbation amplitude the 
latter has value about the period of the harmonic field. 

IV. THE BROADBAND PERTURBATION 

The threshold character of the separatrix crossing in the harmonic field stems from 
the termination of the resonant energy absorption before the value E s is achieved. If the 
perturbation is a broadband noise, then the energy absorption can go on indefinitely at 
arbitrarily small field amplitude. The problem of the energy absorption from the external 
noise has been studied extensively as a part of the theory of the escape over the potential 
barriers in dissipative systems in contact with the heat bath. The averaged evolution of 
the system can be described as a process of diffusion on the energy axis. The equation 



that governs this process could be derived from the Fokker - Planck equation p4| , p!5| . It 
is inconvenient, however, to adjust the known results to our case since our Hamiltonian 
(7,8) has rather unusual structure. Instead we derive the equation for the energy diffusion 
considering the system formally as quantum one, starting from quantum kinetic equations 
and going to the classical limit h — > in the result (cf. [0,01). 



Let's consider a quantum system with the unperturbed Hamiltonian H with one degree 
of freedom and discrete energy spectrum under the perturbation (t) where V depends 
on the dynamical variables of the system and £ (t) is a stationary weak broadband noise 
specified by its spectral density S(u). The state of the system could be described by the 
probabilities p n of finding it in the quantum state \n). The evolution of these probabilities 
obeys the system of master equations 

dp n 00 • 00 

= ~Pn W n,n+k + ^ Pn+kW n+k>n . (21) 

k=—n k=—n 

The rates of transitions W n+k ^ n are determined by the perturbation theory formula 

W^n.n+fc = -^2 |Ki,ra+fc| S (— LO n ^ n+k ) , (22) 

where Ki.n+fc are the matrix elements of the perturbation and S (— w n ,n+fc) is the spectral 
density of noise at the frequency of transition. Let's take the probabilities to be functions 
not on the level number n, but on its energy: p n — > p(E n ). In the quasiclassical case the 
energy spectrum of the system could be related to the frequency of its classical motion at a 
given energy f2 (E): 

E n+k = E n ± htt {e u ± k-rPj , (23) 

and the matrix elements of the perturbation V 1hn+k could be replaced by the Fourier compo- 
nents of the unperturbed motion of the dynamical variable that corresponds to the operator 
V: if 

00 

V(t)= Yl V k e~ ikQt , (24) 

k=— 00 

then 

V n , n+k ^ V k (e u + k^pj . (25) 

We assume p (E) to be a smooth function, and expand its value to the terms of the second 
order in H and the values of E n+k and V njTl+k to the first order in h. After the substitution 
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of these expansions in Eq. (21) and going to the limit h — > we obtain purely classical 
equation. At this point we restrict our consideration to the case of the white noise with the 
constant spectral density S (cu) = 1. For this case we have 



dp = Q _d_ 
dt dE 



(26) 



where the energy diffusion coefficient D (E) is given by the expression 

D (E) = 2-kQ (E) k 2 V% (E) = - j — dt. (27) 

k=l 

The energy dependence of the diffusion coefficient D for the unperturbed system (7) with 
the perturbation V(A, 0) = Vl — A 2 cos(0/2) is shown in Fig. 5. We note a discontinuity 
of D(E) at the separatrix value of energy: D(E S + 0) = 2D(E S — 0). This jump is not a 
direct consequence of the presence of the separatrix, but reflects both the global structure 
of the phase space and the behaviour of the perturbation V in the neighbourhood of the 
separatrix. 

Since the classical distribution in energy w (E) is connected to the probability density 
p (E) by the relation w (E) = p (E) [HQ (E)]' 1 , we have the equation for the w (E) in the 
form 

«H^H- (28) 

The stationary solution of Eq. (28) has the form given by Eq. (18), but the second line 
holds now in all the range E s < E < E + . If the system initially was in one of the stable 
states, then in course of the energy diffusion it has an opportunity to migrate to vicinity 
of the contrary stable state. The characteristic time of the average transfer is determined 
by the relaxation time of the distribution to its stationary value; from Eq. (26) it can be 
estimated as 

(E+ - E„f 

T ~ (QD) ' (29) 
where the angular brackets denote averaging over the energy. 
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V. CONCLUSIONS 



The main question addressed by this paper is: "If the system, that is described by the 
one-dimensional non-linear Scrodinger equation with the potential of the symmetric double 
well, is initially in one of the lowest asymmetric (self-trapped) states, then can the time- 
dependent perturbation transfer the system completely to the opposite asymmetric mode?" 
The answer to it comes to be "Yes, but almost completely and only by chance". 

For the harmonic perturbation with the amplitude F that exceeds the threshold F c (ui), 
the system's motion is captured in the stochastic layer that embraces both domains of li- 
bration around the stable states and a strip around the separatrix. When moving in this 
domain, the system can come arbitrarily close to the opposite asymmetric state. However, 
the nature of this process is purely chaotic and, thence, unpredictable. There is no way 
to create for the nonlinear system the "ir - pulse" that will transfer the system un- 
ambiguously from one of the stable states to another. Lastly we note, that the threshold 
magnitude of the perturbation is only numerically small in comparison with the depth of the 
self-trapping wells: to make the transfer possible, the system must be perturbed strongly. 

For the system under the influence of the white noise (or, generally speaking, any suffi- 
ciently broadband noise) the process of the energy diffusion eventually spreads the probabil- 
ity density over all phase space of the system H . In this case the system can occasionally 
come close to the opposite stable state. We note, however, that the probability of finding 
the system within one of the self-trapping wells is rather small; with our standard set of 
parameters it is only about few percent. 

The main approximation of our calculations consisted in the truncation of the expansion 
of the wave function to just two modes (see Eq. (4)). It was justified by the high quality of 
this approximation in representation the unperturbed self-trapped states 0. Whether this 
accuracy will hold for the seriously perturbed system is quite a different question. For the 
harmonic perturbation of moderate amplitude the system stays locked within the narrow 
energy domain (see Fig. 4), and the addition of contributions of modes $j with i > 2, that 
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will lead to the extension of the energy space of the system, will have little influence. 

The situation may be different for the broadband noise, where the system can reach any 
point in the phase space. However the main contribution to the energy diffusion coefficient 
comes from the frequencies that are lower than fl : in particular, at the separatrix they 
contribute about 0.66 of the total value. Thus if the spectrum of noise has cutoff of high 
frequencies just above the flo, then the energy diffusion ceases at the energy Eh > E s at 
which Q(E h ) = (for our set of parameters E h = —0.356) and the system stays locked 
within the restricted energy domain. Then in parallel to the case of the harmonic field we 
can conclude on the unimportance of the extension of the energy space by addition of higher 
modes. 
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FIGURES CAPTIONS 



FIG.l. Phase trajectories of the system (7) in the absence of perturbations on the plane 
A - for different values of energy: A - E = -3.7, B - E = E s = -3.487, C - E = 0, D - 
E = 5. 

FIG. 2. Dependence of threshold amplitudes F c of the non-diagonal harmonic perturbation 
(10) on the relative frequency co/flo for the nonlinear Bloch equations (7) with the initial 
phase = 0. 

FIG. 3. Dependence of threshold amplitudes F c of the perturbation on the frequency uo for 
the Duffing oscillator (11). The dashed line - the zeroth (linear) approximation, o - estimates 
found from Eqs. (14) (solid line is an eyeguide); • - results of the numerical experiment. 

FIG. 4. Distribution w of system's values of energy under the influence of the over-threshold 
harmonic perturbation with u = Q = 2.887 and F = 2F C = 0.2. Dashed line -theoretical 
distribution Eq. (18), solid line - results of the numerical experiment. 

FIG. 5. Dependence of the coefficient of the energy diffusion D under the influence of the 
white noise of the unit spectral density on the energy E. 
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